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Abstract 

The nucleon electromagnetic form factors are calculated in light 
cone QCD sum rules framework using the most general form of the 
nucleon interpolating current. Using two forms of the distribution 
amplitudes (DA's), predictions for the form factors are presented and 
compared with existing experimental data. It is shown that our results 
describe remarkably well the existing experimental data. 
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1 Introduction 



The nucleoli electromagnetic (EM) form factors are the fundamental ob- 
jects for understanding their internal structure. The internal structures of 
the nucleon are usually described in terms of the electromagnetic Dirac and 
Pauli form factors Fi{q^) and -^2(0'^) or equivalently the electric and mag- 
netic dipole Sachs form factors Ge{(i'^) and Gm{,(1^), respectively (for a recent 
status of experiments and phenomenology of the form factors see |T]). 

Until a few years ago, the nucleon electromagnetic form factors are studied 
in unpolarized elastic electron-nucleon scattering through a virtual photon 
exchange. It is shown in the pioneering work [2J that the polarization effects, 
i.e., scattering of polarized electrons from polarized target, can play essential 
role for a more accurate determination of the nucleon electromagnetic form 
factors. The main result of [2] is that, unlike the unpolarized elastic cross 
section, which is proportional to the sum of squares of the form factors, the 
polarized cross section contains also interference terms of the form factors 
Ge{ci'^) and GM{<f)- Studying various polarization observables allows more 
accurate determination of these form factors. 

Recent developments in experimental instruments allow to produce polar- 
ized electron beams and polarized protons, which gives the opportunity for a 
more precise separation of the GE{<f) and GM{q'^) form factors. The electron- 
proton scattering experiments, which are performed at Jefferson Laboratory 
using the polarized electrons and polarized proton, show strong deviation 
from the theoretical predictions [Sj HI [5l [6] , i.e., the ratio F2(g^)/Fi(g^) does 
not behave as is expected from previous experiments and as is predicted by 
the perturbative QCD (for review see [7] and references therein) . For under- 
standing this unexpected result, some model- independent non perturbative 
method is needed. Among all existing nonperturbative approaches QCD sum 
rule is more attractive and powerful, because it is based on the fundamental 
QCD Lagrangian. 

The goal of our work is the calculation of the electromagnetic form factors 
of nucleon using the light cone QCD sum rule (LCQSR) and most general 
form of the interpolating current for nucleon. In this approach the form 
factors of the nucleons are expressed in terms of distribution amplitude of 
the nucleon. Note that, this problem is investigated for the loffe current in 
the framework of the LCQSR in [8] and the traditional sum rules in [9]. In 
|10j . an improved version of the ChernyakZhitnitsky current is used. The 
paper is organized in following way. in section 2 , we present the result for 
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the nucleoli electromagnetic form factors in the LCQSR method. Section 3 
is devoted to the numerical analysis, discussion and conclusion. 



2 Electromagnetic form factors of nucleon in 
LCQSR 

In this section EM form factors of miclcon are calculated within the light cone 
QCD sum rules method. The electromagnetic form factors of nucleon are 
defined by the matrix element of the electromagnetic current between the 
initial and final nucleon states {N{p') \ \ N{p)). The most general form 
of this matrix element satisfying the Lorentz invariance and electromagnetic 
current conservation is 



(iV(p') I Jt{Q) I N{p)) = N{p') 



ZniN 



N{p), 

(1) 



where — —q^, is the negative of the square of the virtual photon mo- 
mentum, q — p — p' and Fi and F2 are the Dirac and Pauli form factors, 
respectively. 

Another set of nucleon form factors is the so called Sachs form factors, 
which are defined in terms of the Fi{Q'^) and -F2(Q^) as follows: 

G'M(Q') = Fi(g2) + F2(g2), 



Ge{Q-')^F,{Q-')-^^F,{Q-% (2) 



At the static hmit, values at = q are G^(0) = 1, G^(0) = 0, G^(0) = 
= 2.792847337(29) and ^^^(0) = = -1.91304272(45), where fip and 
Hn are the anomalous magnetic moments of the proton and neutron in units 
of the Bohr magneton. 

After these preliminary remarks, we proceed to calculate the electromag- 
netic form factors of nucleon in LCQSR. The basic object of the LCQSR is 
a suitably chosen correlation function. In this study, it is chosen as: 

n,(p,g) = t J d'xe'^^iO I T{J^(0)J^'(a;)} | N{p)), (3) 
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which describes the transition of the nucleon N{p) to the nucleon N{p — q) 
via the EM current. The interpolating current for the nucleon is chosen as 



where A\ = I, A\ = A\ = 75, = P, and C is the charge conjugation 
operator, and a, b, c are the color indices. The electromagnetic current is: 

J^\x) = CuU'jxu + edd'jxd, (5) 

and the choice /3 = — 1 corresponds to the loffe current. The main idea of the 
LCQSR method is to calculate the correlation function in terms of the form 
factors at hadron level, as well as in terms of the quark and gluon degrees 
of freedom. Equating two representations of the correlation function and 
performing a Borel transformation in order to suppress the contributions of 
the higher states and continuum, we get sum rules for the EM form factors 
of the nucleon. 

Let us first calculate the physical part of the correlator Q. The contri- 
bution of the nucleon to the correlation function is given by 

n rp,^) - y <Q I '^'^(Q) I N{p',s)){Nip',s) I (0) I Njp)) ^ 
g m% — p'"^ 

The matrix element (0 | J^(0) | N{p' , s)) in ([6]) is determined in the following 
way: 

(0| J^(0)|iV(p',.)) = A^iV(p',.), (7) 

where Aa? is the coupling constant of the nucleon to the current J^(0). The 
matrix element {N{p',s) \ J|^(0) | N{p)) is parameterized in terms of the 
form factors Fi and F2 via Eq. ([T]). Summing over spins of the nucleons 

Y,N{p',s)N{p',s) =/ + m^, (8) 

s 

and using Eqs. ([T]), (Ej) and (JTj), we obtain the following expression for the 
contribution of nucleon to the correlation function 



nA(p, q) = A i'' + 



N{p) + ■■■ 

(9) 
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where • • ■ stand for the contributions to the correlation functions from the 
higher states and continuum. It follows from expression that, the cor- 
relation function contains numerous structures and in principle all of them 
can be used in determination of the electromagnetic form factors of nucleons. 
In further analysis, we choose the independent structures containing p\, and 
Pa^ for obtaining Fi and F2, respectively. 

The theoretical part of the correlator can be calculated in LCQSR in 
deep Euclidean region p'"^ = {p — « in terms of the nucleon DA's. 
These nucleon DA's for all three quarks have been studied in great detail in 
[H [m [To] . Using the explicit expression for the currents and carrying out all 
contractions, the correlation function takes the form 

(Ha), = iJdV^p^l^ 

euiCAiU[AiSui-xhx\ 4e"''^(0|<(0)4(x)rf^;(0)|Ar(p)) 



+ euiAl),^[iCA\y Sui-xhx\Je'^'^OK{0)ul{x)d;{Q)\N{p)) 

(10) 

in X representation, where A is a Lorentz index, and a, 7, p and are spinor 
indices. S{x) is the light cone expanded light quark full propagator [12] 
having the form: 

i ^ mgX^ . f 7 r ^ 





- ^^'^-^^^ii^]' (11) 

where ml = (0.8 ± 0.2) GeV^ and G^j, is the gluon field strength tensor. 
The terms proportional to the gluon strength tensor can give contribution 
to four- and five-particle distribution functions but they are expected to be 
very small [H Ull [10] and for this reason we will neglect these amplitudes 
in further analysis. The terms proportional to < gg > can also be omitted 
because Borel transformation eliminates these terms and hence only the first 
term in Eq. ( |TT|) is relevant for our discussion. It follows from Eq. ( [TOi) that 
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for the calculation of Ilx{p, q) we need to know the matrix element 



(0 I Ae'''%l{aix)ul{a2x)d;{asx) \ N{p)). (12) 

It is shown in [11] that the general Lorentz decomposition of this matrix 
element is symmetric with respect to interchange of the momentum fractions 
of the u-quarks: 

(0 I 4e'^''X(«ia;)<(a2a:)rf^(a3x) | Nip)) = ^ KVf iT2N{p)y , (13) 

where N{p) on the right is the nucleon spinor, ri^2 are certain Dirac structures 
over which the sum is carried out, are positive numbers which satisfy 
ai + 02 + as = 1, and K are the distribution amplitudes, depending on eight 
nonperturbative parameters. Explicit expressions of all DA'S and the values 
of eight nonperturbative parameters can be found in [S], [TTl [TUl [13] • 

Omitting the details of calculations of the theoretical part, choosing the 
coefficients of the structures px, and px ^, equating both representation of 
the correlation function and applying the Borel transformation with respect 
to the variable p'^ = {p — q)"^, which suppress the contributions of the higher 
states and continuum, we obtain following sum rules for the form factors Fi 
and F2: 



2Xj\i I J to Jo 



2n,,Mx.^{i-P) 



+4(7^l7(x^)-27^l9(x,))(l + /3) 



1 — X2 



+e^mr, / dx2 I dx, H ^e'^^'^'^')/^^^ 

to "'0 Jto ti 



Ml 



H2o,-i8(a;.)(l + /3) - n^{xi){-l + (3) 

{2H2oAxi){l + I3){Q^ + s(ti, Q') + ml{-l + t^))} 



{Hl5,-14(Xi)tl(l - /3) - 47i:21,24(Xi)tl(l + (3) 
+2nio{xi){-l + /?)(ti - X2) + 2(7^i6(Xi)(-l + 13) + 27^24(X.)(1 + P))X2} 



euUiN / dx2 

Jto JO 



1 — X2 



dxie-"'/^^B. 



27^2o,l8(a;^)(l + /?)(Q' + S0 + m%i-l + to)) 



+ m%tl 

{7^_8,9(X,)(1 -13)- (37^21,24 (X.) + 87^23(X.))(1 + P)} to 

+2Hio(xi)(-l + P)ito - X2) + 2{nie{xi){-l + P)+ n2i{xi){l + P))x2 



+mlf 



+edv'^{Q^^3) + eMQ^P)], 



X 



N 



{cu I dX2 I 
Jto Jo 



X2 



-e„m^ / dx2 / dxi 

Jto Jo Jti 



^2 dtl 



-s(ti,Q2)/M| 



to ti 



Ml 



nsMxi){l -13) + 2(Hi8,2o(a;0 + 2n2iM^i) + m23ixi)){l + p) 



n22{Xi){l + P)x2 
1—X2 



dxo 



dx^e-'^'^B 



to JO \Q'^ + m%tl 

-2{nisA^i) + 2:^21,22(^0 + 4H23(^i))(l + l3)tQ + 47^22(2:0(1 + P)x2 
+edV2{Q\ P) + euV2{Q\ P) 
where 



J^{Xi) = J^{xi, X2, 1 - Xi- X2), 
T{x[) = ^(Xi,l -Xi -X3,X3), 

y 

with ^0(^0, Q^) being solution of the equation s{to, Q^) — Sq, and 



(16) 



r)i{Q'^,P) = niN \ / dx3 



to 



l-X-3 



(^0 



(15) 



(Hl,l7,3(^0-2Hl9(x,))(l+/9) 



+Hi3,7(a;0(-1 + /?) 

1 



+ 



Jto Jo Jto 



M% 



to JO Jto 
2 / _2 I r^'i I /^2 



-n22{xi)m%{-m% + + s(ii, g'))(l + /3)x3 

H22(2:i)"^^("^^(-i + 2ti - 2x3) + + s(ti, g'))(i + p) 



2Mlti 



(m 



N 



- sih, g2)) {{Hisixi) - 37f2o(a;^))(l + P) + 2H^Ax,){-l + 



1 



Ml 



rUj^ {^-12,15 -6,9(2;i)(l - /?) + (^18 -2,24,4,21 (a^i) + 27i_20,-22,23(2;i)) (1 + /?)} 



1 



^i2,6(a;.)(l - /3) + 7f-i8,2o(a;i)(l + /?) 
1 



1 



Jto Jo 

n22{xi)m%{l + P)tl{Q^ + S0 + m%{-l + 2to)){to - Xs) 



1 



{Q' + m%tir 

+ m^(2to(-l + 3to - 2x3) + X3)) 
1 



2n22(xi)m%(l + p)4{Q' + S0 + m^(-l + 2to))(to - X3) 
n22{xi)ml{l + (3)tl{{Q^ + So)(2to - X3) 



+ 



1 



2{Q^ + m%tl) 



rl rl—xs 

V2(Q^P)= / dxs / dxie-^(^3'« 

Jto Jo 



2)/M| 



X3 



+mN 



Jto Jo Jti 



'to 



n22{xi)m%{l + p) 



+ 



1 



3 



H22{xi)mN{l + - + ^^(1 + 2x3)) 

n22{xi)mN{l + /3)(g^ + s(ti, Q^) - ml)x3 



'H22{Xi)mN{l + (3)X3 



+ (^2,-20,-21,-4(2:0 + m22{Xi) - 27i23(a;i))(l + 13) 



2n22{xi)ml{l + P)tlx3 

-H2o(a:.)(l + /3)to{3(g' + so) + m^(-3 + 4to)} + 27f6,i2(a;i)(-l + /9)(g' + so 
+m^(-l + to))to + 2n24{xi)m%{l + f3)to{to - x^) 
+2nisixi)il + PMQ^ + S0 + m%i-l + 2to)) 

+2mlHg,-i5{xi){-l + (3)tl + 2m^(7f4,-2,2i(a;i) + 2n23,-22{xi)){l + p)tl 



(17) 
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+mN 



Jtn Jo 



rriN 



+ 



+ 



MUQ' + m%tl 
1 



(Q2 + ^^^2)3 
1 



(Q2 + ^2^^2)2 

rriN 



n22{xi){l + P){Q^ + so + m%{-l + 2to)){to - xs) 
-2H22{xi)m%{l + p)tl{Q^ + So + m^(-l + 2to)){to - x^) 
n22{xi)mN{l + P)tl{Q^ + So + ml{-l + 4to - 2x3)) 

-12,-15,-6,9(3:^1) (1 — P) + ("^2,-20,-21,22,-4 (2^0 — 27^23(3:^1)) (1 + P)to 



H22 (2^0 (1 + /9) (3^:3-2^0 



and r)'^{Q'^, /3), (i = 1, 2) are obtained from r}i{Q^, (3) by replacing 0:3 with X2 
and replacing T{xi) with J-'{x'j) in the integrals. In the above equations, we 
have used the short hand notations for the functions Ti-±i,±j,... = ±?ij ±7^^..., 
and Tii are defined in terms of the distribution amplitudes as follows: 



Hi — Si 


7^2 — 'S'1,-2 


Hz^Pi 


7^4 = Pl,-2 


n.^Vi 


^6 — ^1,-2,-3 




7^8 = -21/1,-5 + 14,4 


H9 — 


7^10 = —14,-2,-3,-4,-5,6 


Hn = Ai 


7il2 — — ^1,-2,3 


Tiia — A3 


7^14 = — 2Ai,_5 — ^43,4 


7^15 = ^3,-4 


"^16 — ^1, -2,3,4,-5,6 


7^17 = Ti 


7^18 = 7i,2 — 2T3 


'Hig — T-j 


7^20 — Ti-2 — 2T-J 


7^21 = ~21__5 + 2^8 


TC22 — 72,-3,-4,5,7,8 


7^23 = 27,-8 


7^24 = —7l, -2,-5,6 + 2T7, 



(19) 

where for any distribution amplitudes, X-tj -tj, = ±Xj ±Xj... are also used. 
The overlap amplitude of the nucleon interpolating current with nucleon is 



(18) 
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determined from sum rule and its expression is 
A 



^IIMI j^E,{x){5 + 2(3 + P') + 



2 ^ ml/Ml 



\2567r4^ 



6 



- 6(1 - p^) <dd> 



+ (-1 + <uu> 
+ (-1 + pf <uu> 



rrin 



24M| 



< uu > 



12(1 -Z^^) <dd> 



(20) 



where x = sq/M"^ and the function 



En{x) = l-e-^±^ (21) 

k=l 

corresponds to the continuum subtraction. 

3 Numerical results 

It follows from explicit expressions of the sum rules for the nucleon electro- 
magnetic form factors that, the nucleon DA's are the principal input parame- 
ters, whose explicit expressions can be found in fS]. These DA's contain non- 
perturbative parameters which should be determined in some framework. In 
the present work, we consider two different determination of these input pa- 
rameters: a) All eight nonperturbative parameters /at, Ai, A2, Vf, A", ff, /" 
and /2 are estimated within QCD sum rules method [8],|TT1|10] (setl), b) The 
condition that the next to leading conformal spin contributions vanish, fixes 
five of the eight parameters. This is the so called asymptotic set. The values 
of all nonperturbative parameters are (see [S]): 



fN 

Ai 
A2 



= (5.0 ± 0.5) X 10"^ G'e1/^ 
= -(2.7 ± 0.9) X 10-2 GeV^2 
= (5.4 ± 1.9) X 10^2 G'e1/^ 



(22) 
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setl 



asymptotic 
^1=0, 



A 



u 



1 



0.38 ±0.15, 




■u 



0.07 ±0.05, 



0.22 ±0.05, 



0.23 ±0.03, 



0.40 ±0.05, 




(23) 



The continuum threshold that appears in the continuum subtraction is 
determined from the mass sum rules as sq = 2.25 GeV^. There are two auxil- 
iary parameters of the sum rules: the Borel parameter and the parameter 
f3. The Borel mass square is the artificial parameter of the sum rules 
and therefore wee need to find a region of M^, where physically measurable 
quantities, in our case electromagnetic form factors, be independent of M^. 
Lower bound of is determined from condition that contribution from 
higher states and continuum in the correlator should be enough small, upper 
bound of is determined from condition that series of the light cone ex- 
pansion with increasing twist should be convergent. Our numerical analysis 
shows that both conditions are satisfied in the region IGeV"^ < < 2GeV^, 
which we will use in numerical analysis. 

The other auxiliary parameter f3 is chosen in a region such that, the 
predictions are independent of the precise value of (3 in that region. In our 
analysis, it is shown that in the region —0.5 < cosO < 0.5 the form factors are 
practically insensitive to the variation of (3, where 6 is defined as tan6 = 13. 
Note that the analysis of mass sum rules and magnetic moments of octet 
baryons [11] leads to the very close region for cosO, i.e. —0.6 < cosO < 0.3. 
Also, it is observed in [15] that the optimal value of /5 is /3 = — 1.2(cos^^ = 
—0.64), which follows from the Monte Carlo analysis. 

In Fig. [1], we present the dependence of the proton magnetic form factor 
GIj/hpGd on at sq = 2.25 GeV^, M| = 1.2 GeV^ for two sets of DA's, at 
fixed values of parameter f3. In this figure, we also present the experimental 
results [ISl [13 dB]. From this figure, we see that the dependencies, as 
well as the magnitude of proton magnetic form factor are rather in good 
agreement with the experimental data, especially for the set 1 of DA's and 
loffe current {f3 = —1). The dependence of the ratio of the proton electric 



10 



form factor to the magnetic form factor ^pG^^/G\j on at Sq = 2.25 GeV^, 
= 1.2 GeV^ for two sets of DA's, at fixed values of parameter (3 is 
depicted in Fig. [2l From this figure it follows that, practically, both sets of 
DA's well describe the existing experimental results, except for /? = 5 and 
/3 = —1 of set 1. For large values of Q^, > 4 GeV"^, the experimental 
results obtained in [1] and in [T7l[IB] are not in agreement. Whereas (3 = —1 
describes better the data in larger values of describe better the data 
in [H]. 

The LCQSR results for the neutron magnetic (normalized to the dipole 
form factor ) and electric form factors are given in Fig. [3] and Fig. HJ 
respectively. From Fig. [21 we see that the magnetic form factor of neutron 
reproduce experimental data very well at /3 = — 1 for both sets of DA's. 
Neutron electric form factor is in a good agreement with the experimental 
result for all cases. 

Analysis of the experimental results (for review see [7] and references 
therein) lead that the magnetic form factors of the nucleon are very well 
described by the dipole formula 

^jvf (Q^) = 7 ^~r^ 72 = f^n,pGD- (24) 



(1 



(0.71 GeVy 



The measured values of the electric form factors of the neutron are given in 
[201 [21]. 

In [22l [23] , the following large behavior of the electromagnetic form 
factors is obtained 

where A = 300 MeV. In Fig. [5] ([6]), we present the logarithmic scale pre- 
diction, i.e. (l/15)/n-2(QVA2)Q2^2(Q^)/i^i(Q^) for the proton (neutron), 
with available experimental data [24J at fixed values of (3 for two sets of 
DA's. From these figures, we see that our prediction for the proton for 
In^^^Q"^ / A'^)Q'^F2{Q^)/ Fi{Q^) is in good agreement with experimental data 
except for /3 = — 1 case for both DA's, and /3 = —5 case for setl. For the 
neutron case only setl for /? = — 1 describes quite successfully the existing 
experimental data. 

Finally, in Fig. [TJ as an example on the dependence of the predictions on 
/3, we present the dependence of proton magnetic form factor normalized to 
the dipole form factor G\^/ UpGr, on cosO, for both sets of DA's at two fixed 
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values of Q^. It follows from this graph that, in the chosen region of /3, i.e. in 
the region —0.5 < cosO < 0.5, the form factor G\j is practically insensitive 
to the variation of j3. 

In conclusion, in present work, we calculate the nucleon electromagnetic 
form factors using the most general form of the nucleon interpolating current 
in the light cone QCD sum rules. The sum rules for these form factors are 
obtained. Using two forms of the DA's, we calculate sum rules predictions for 
these form factors and compare them with existing experimental data. We 
obtain that our results are in a good agreement with the existing experimental 
data. More precisely, at different values of /3, our results for the form factors 
reproduce the experimental data. Finally, we obtained the "working region 
for 

Our final remark is that in order to answer to the question which (3 is 
more preferable, both theoretical and experimental studies have to be refined. 
Prom theoretical part 0{as) corrections to the distributions amplitudes and 
more accurate determination of the DA's are needed. Prom experimental 
data, the discrepencies between various data has to be eliminated. 
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Figure 1: The dependence of G^j/fipGo on at sq = 2.25 GeV^,M^ = 
1.2 GeV^ for f3 = —1, — 5 and 5. The boxes correspond to experimental 
data in [16], the diamonds to [17] and the up-triangles to [18] . The lines 
with circles correspond to setl and the hnes without any circles correspond 
to the asymptotic DA's 
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Figure 2: The same as Fig. [H but for upG^jG^. The boxes/diamonds/up- 
triangles/down-triangles/right-triangles/left-triangles correspond to experi- 
mental data given in [16]/[T7j/[6]/[l]/[5]/[3] respectively 
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Figure 3: The same as Fig. [T]but for G\j j jinGo- The boxes correspond to 
experimental data ([19]) 
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Figure 4: The same as Fig. [T] but for G^. The boxes are correspond to 
experimental data ([19]) 
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Figure 5: The same as Fig. [Jbut for Q^ln-'^{^)F^ / Ff where A = SOOMeV. 
The boxes correspond to experimental data ( [24\ ] 
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Figure 6: The dependence of j^QHn-^{^)F^ / on at sq = 
2.25 GeV^,M]^ = 1.2 GeV^,A = SOOMeV. The boxes correspond to ex- 
perimental data ([21]) 
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Figure 7: The dependence of G^/ixpGo on cos9 at Sq — 2.25 GeV'^.M'^ = 
1.2 GeV^ for two different values of i.e. = 2 GeV^ and = 4 Ge^^ 
The hues with circles correspond to setl and the lines without any circles 
correspond to the asymptotic wavefunctions 
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